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We propose a technique to measure multi-spin correlation functions of arbitrary range as deter- 
mined by the ground states of spinful cold atoms in optical lattices. We show that an observation 
of the atomic version of the Stokes parameters, using focused lasers and microwave pulsing, can be 
related to n-spin correlators. We discuss the possibility of detecting not only ground state static 
spin correlations, but also time-dependent spin wave dynamics as a demonstrative example using 
our proposed technique. 



PACS numbers: 03.75.Lm, 75.10.Pq, 03.75.Mn, 39.25.+k 



I. INTRODUCTION 



The advent of optical lattice confinement of ultracold 
atomic gases P, d, d, Q opens the possibility of observ- 
ing a vast array of phenomena in quantum condensed 
systems Q. In particular, optical lattice systems may 
turn out to be the ideal tools for the analog simulation 
of various strongly correlated interacting lattice models 
(e.g. Hubbard model [1, Q, Kitaev model |6]) studied in 
solid state physics. The great advantage of optical lat- 
tices as analog simulators of strongly correlated Hamil- 
tonians lies in the ability of optical lattices to accurately 
implement lattice models without impurities, defects, lat- 
tice phonons and other complications which can obscure 
the observation of quantum degenerate phenomena in the 
solid state. 

In this context optical lattices can support a variety 
of interacting spin models which to date have been only 
approximately or indirectly observed in nature or remain 
as rather deep but unobserved mathematical constructs. 
Three exciting possibilities are currently the subject of 
active study 0. The first (and the most direct) envi- 
sions simulation of conventional condensed matter spin 
lattice models in optical lattices. Quantum magnetism 
arising from strong correlation leads to many-body spin 
ground states that can be characterized by spin order pa- 
rameters. Spin order can, in some cases, show long range 
behavior arising from spontaneous symmetry breaking, 
e.g. ferromagnetism and antiferromagnetism. Such long 
range spin ordering phenomena are reasonably well un- 
derstood in most cases. Recent work also relates conven- 
tional spin order parameters to entanglement measures 
which yield scaling behavior near quantum phase transi- 
tions [3, Q • The second possibility, simulation of topo- 
logical spin states, arises from the surprising fact that 
optical lattices can also (at least in principle) host more 
complicated spin models previously thought to be aca- 
demic. The ground states of these models do not fall 
within the conventional Landau paradigm, i.e. there is 
no spontaneously broken symmetry, but show topological 
ordering and, as a result, display nontrivial short range 
behavior in spin correlation functions. Examples include 
the chiral spin liquid model and the Kitaev model 



[1, [T^, [Tl|, [T2| ■ And finally, optical lattices are also par- 
ticularly well suited to realize coherent and collective spin 
dynamics because dissipation can be kept to suitably low 
levels [il. 

While optical lattices offer the possibility of realizing 
all three of the above examples one glaring question re- 
mains. Once a suitable spin Hamiltonian is realized, how 
do we observe the vast array of predicted phenomena in 
spin-optical lattices? To date time of fiight measurements 
have proven to yield detailed information related to two 
types of important correlation functions of many-body 
ground states of particles trapped in optical lattices. The 
first, a first order correlation function (the momentum 
distribution) . indicates ordering in one-point correlation 
functions The second is a second order correlation 
function (the noise distribution) which indicates ordering 
in two-point density-density correlators [l^, [13, [3 ■ 
The former can, for example, detect long range phase 
coherence while, as we will see below, the latter is best 
suited to probe long range order in two-point correlation 
functions, e.g. the lattice spin-spin correlation function. 
We note that recent proposals suggest that time of flight 
imaging can in principle be used to extract other corre- 
lation functions [l^, [23| . 

In this paper we propose a technique to observe equal 
time n-spin correlation functions characterizing both long 
and short range spin ordering useful in studying all three 
classes of spin lattice phenomena mentioned above. Our 
proposal utilizes realistic experimental techniques involv- 
ing focused lasers, microwave pulsing and fluorescence 
detection to effectively measure a general n-spin correla- 
tion function deflned by: 



^}^(vI/|I]a;- I*) 



(1) 



k=l 



where ^' is the many-body wavefunction of the atomic en- 
semble, {jk} is a set of sites, and cr^^'' {aj^ = 0, 1, 2, 3) are 
Pauli spin operators at sites jk with the notation a'-* = I, 
= (7^, = (7^, and cr" = cr^. Examples of order 
detectable with one, two and three-spin correlation func- 
tions are magnetization ((o"!) — 1), anti-ferromagnetic 

order {{'^jO'j,') — (— l)-'^-' ), and chiral spin liquid order 
(((Tj • (cTj' X <7j")) — 1), to name a few. 



In general our proposed technique can be used to ex- 
perimentally characterize a broad class of spin-lattice 
models of the form: 



H{J-A) = J{t)Y, 



A 



M 
k=l 



where J has dimensions of energy and can vary adiabat- 
ically as a function of time, while the dimensionlesf 
parameters ^{jj.} are kept fixed. For example, M = 2 
represents the usual two-body Heisenberg model. Severa 
proposals now exist for simulating two-body Heisenberg 
models In the following we, as an example, con- 

sider optical lattice implementations of the Heisenberg 
XXZ model: 



E 



(3: 

where (j, j') denotes nearest neighbors and A and J are 
model parameters that can be adjusted by, for example, 
varying the intensity of lattice laser beams [lQ| . 

The paper is organized as follows: in Section II we 
show that, in practice, short range spin correlations are 
difficult to detect in noise correlation measurements. Sec- 
tion HI lays out a general procedure for detecting n-spin 
correlations with local probes. An experimental scheme 
for realizing such a procedure is proposed and a quan- 
titative feasibility analysis is presented. Section IV is 
devoted to an investigation of the time-dependence of 
correlation functions using the technique. In particular, 
we show how the technique can be used to engineer and 
probe spin wave dynamics. We conclude in Section V. 



II. LOCAL CORRELATIONS IN TIME OF 
FLIGHT: 

We first discuss the measurement of spin-spin cor- 
relation functions by analyzing spatial noise corre- 
lations (two point density-density correlations, i.e. 
((n (r) n {f')) — {n {r)) {n {r'))) / {n (f)) {n {r')) in time of 
flight images from atoms confined to an optical lattice 
modeled by the XXZ Hamiltonian. The ground states 
of this and a variety of spin models can be characterized 
by the spin-spin correlation function between different 
sites. For instance, the spin-spin correlation function in 
a one dimensional XXZ spin chain (with J > 0), shows 
power-law decay 



(a|4)^(-ir^ /\j-f\ 



(4) 



in the critical regime (—1 < A < 1), where 77 = 
l/(l — icos^^A). In principle this correlation func- 
tion can be probed by spatial noise correlation in time of 
flight. 

We argue that, in practice, short range correlations 
(e.g. 77 > 1 in the XXZ model) are difficult to detect 
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Figure 1: (Color online) Noise correlation plotted as a func- 
tion of wavevector of the one- dimensional XXZ model. The 
solid (dashed) line corresponds to a ground state with 77 = 
0, long-range {jj — 2, short-range) spin correlator. The ampli- 
tudes are normalized by the maximum for anti-ferromagnetic 
order {rj — 0). The number of atoms in panel a (b) is A'' = 20 
{N = 200). 



in time of flight noise correlation measurements. To see 
this note that the noise correlation signal is proportional 
to fisj: 



G(Q (r-/)) = 



(5) 



where Q is the lattice wave vector which gets mapped 
into coordinates r and r' in time of flight on the de- 
tection screen, and a is the lattice spacing. Including 
normalization the noise correlation signal is proportional 
to for systems with long range order (e.g. anti- 

ferromagnetic order giving 77 = in the above XXZ 
model) but shows a much weaker scaling for short range 
correlations. In fact the ratio between correlators in a 
ground state with 7; > 1 (short range power law or- 
der) and 77 = (long range anti-ferromagnetic order) 
scales as N^^ making the state with power law correla- 
tions relatively difficult to detect in large systems. To 
illustrate this we compare the calculated noise correla- 
tion amplitude, G, in Fig. [1] for two cases 77 = (solid 
line) and 77 = 2 (dashed line) with A*" = 20 (Fig. la) and 
A^ = 200 (Fig.lb) for the ID XXZ model. We see that 
the correlation amplitude for short range (power-law) or- 
der is extremely small in comparison to long range anti- 
ferromagnetic order for large N . 

The small correlation signal originates from the fact 
that the noise correlation method is in practice a condi- 
tional probability measuring collective properties of the 
whole system, while short range spin correlations de- 
scribe local properties and are therefore best detected 



3 



via local operations. In the following we propose a local 
probe technique to measure local correlations thus pro- 
viding an experimental scheme which compliments the 
time of flight— noise correlation technique, best suited for 
detecting long range order. 



III. DETECTING n-SPIN CORRELATION 
WITH LOCAL PROBES 

A. General procedure 

We find that general n-spin correlators, 
(,{aj^,k ~ I, ...,n}, can be related to the Stokes 
parameters broadly defined in terms of the local reduced 
density matrix p = Tr^^^^ ^.^^^ „} j^*) (vj/l on sites 
{jk,k= l,...,n}, where the trace is taken on all sites 
except the set {jk}- The Stokes parameters [2l| for 

the density matrix p are Saj^....aj^ — Tr (^pYYk=i'^'j^ 
leading to the decomposition 



E 



n- 

k=l 



3 k 



(6) 



Using the theory of quantum state tomography [2l|, we 
find the n-spin correlators 

n 

^ {a,, , A: = 1, n} = J] { h".. ) } ± ^ {!<.)}) ' 
fe=i 

(7) 

where the plus (minus) sign indicates a (non-zero) in- 
dex and 1 1 4>aj^ ^ , 0aj^ ^ I denote the measurement ba- 
sis for the atom at jk- We define the measurement basis 
to be: = (|i) + |t))/^/2, \c^i) = (U)-|t))/V2, 

\h) - (li) + ^ |T» /V2, = (|i> - z |T» /V2, 103) = 

li), = It). Finally, P ^} is the probability of 

finding an atom in the state 



The expansion of the product defining ^ then yields a 
quantity central to our proposal: 



^{a,,,fc = l,...,7i} = ^(-l)'p,. 



(8) 



1=1 



where P; is the probability of finding I sites in the states 
|(^^) and n — I sites in \(t)j^)- Eq. ([5]) shows that the 
n-spin correlation function can be written in terms of ex- 
perimental observables. We can now write a specific ex- 
ample of the two-spin correlation function (discussed in 
the previous section) in terms of observables: ^ {3, 3} = 



following subsection we discuss a specific experimental 
procedure designed to extract precisely this quantity us- 
ing local probes of cold atoms confined to optical lattices. 
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B. Proposed Experimental Realization 

We now describe and critically analyze an experimen- 
tal procedure designed to find the probabilities, P;, from 
which we can determine the spin correlation function 



He 



-Jk I 



fc = 1, 



i} through Eq. 



To illustrate our 



technique we consider a specific experimental system: 
®^Rb atoms confined on a single two dimensional {xy 
plane) optical lattice with two hyperfine ground states 
li) = \F=l,mF = -l) and |t) = |P = 2,mF = -2) 
chosen as the spin of each atom. Here the atomic dy- 
namics in the z direction are frozen out by high frequency 
optical traps However, the scheme can be directly 
applied to three dimensional optical lattices by using one 
additional focused laser which propagates along the xy 
plane and plays the same role as the focused laser (prop- 
agating along the z axis) discussed in the following step 
(II). 

In the Mott insulator regime with one atom per lat- 
tice site spin Hamiltonians, defined in Eq. may be 
implemented using spin-dependent lattice potentials in 
the super-exchange limit The spin coupling J {t), 

i.e., the overall prefactor in Eq. can be controlled 
by varying the lattice depth. Our proposed experimental 
procedure will build on such spin systems, although it can 
be generalized to other implementations where H is gen- 
erated by other means. In the spin systems, the ground 
states are strongly correlated many-body spin states and 
their properties can be characterized by the spin correla- 
tions between atoms at different lattice sites. 

Step (I): We start with a many-body spin state $o 
and turn off the spin-spin interactions generated by su- 
per exchange between lattice sites. We achieve this by 
ramping up the spin-dependent lattice depth to ^ bOEn 
adiabatically with respect to the band splitting, where 
En = h? /2m\^ is the photon recoil energy and A is the 
wavelength of the optical lattice. The ramp up time for 
each lattice is chosen properly so that only the overall 
energy scale J {t) in the spin Hamiltonian ([2|) is modi- 
fied, which preserves the highly correlated spin state $0. 
In the deep lattice, the time scale for the spin-spin in- 
teractions (~ h/ J > 10s) becomes much longer than the 
time 1ms) taken to perform the spin correlation mea- 
surement. The spin-spin interactions play no role in the 
measurement process and the following detection steps 
are quickly performed on this "frozen" many-body spin 
state $o- 

Step (II): In this step, we selectively transfer tar- 
get atoms A at site(s) jk (chosen aprior) to a suit- 
able measurement basis {\4>jk) , I't'f^}} from initial states 

;IT)i f' without affecting non-target atoms B at 

other sites (Fig. 2a). Selective manipulation of quantum 
states of single atoms in optical lattices is currently an 
outstanding challenge for investigating physics in optical 
lattice, because spatial periods of typical optical lattices 
are shorter than the optical resolution. In Ref. [23l |. 
a scheme for single atom manipulation using microwave 
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Figure 2: (Color online) Schematic plot of the experimental 
procedure used to measure n-spin correlation functions. A 
indicates one of the target atoms. B indicates all other atoms, 
(a) Target atoms A are transferred to a suitable measurement 
basis using a combination of focused lasers and microwave 
pulses (see step (II)). (b) Atoms at state are transferred 
to state |2) by two microwave 7r-pulses, then atoms at state ||) 
is transferred to state ||) by another microwave 7r-pulse (see 
step (III), (c) Target atoms A axe transferred back to state 
II) from II) (see step (IV)). (d) A detection laser is applied to 
detect the probability of finding target atoms at ||) (see step 
(IV)). 

pulses and focused lasers is proposed and analyzed in de- 
tail. Here we apply the scheme to selectively transfer 
atoms between different measurement bases. 

To manipulate a target atom we adiabatically turn 
on a focused laser that propagates along the z axis having 
the maximal intensity located at A. The spatially varying 
laser intensity / (r) induces position-dependent energy 
shifts 

AS, (r) = A/ (r) (9) 

for two spin states ||) and ||), where the parameter j3i 
for state \i) {i —[ or |) is determined by the focused 
laser parameters. Different polarizations and detunings 
of the focused laser lead to different Pi and thus yield 
different shifts of the hyperfine splittings |Ai?(r)| = 
|Ai?f (r) — Ai?j (r)| between two spin states. We choose 
a cr"*" -polarized laser that drives the 55* — s- 6P tran- 
sition to obtain a small diffraction limit. The wave- 
length \f ~ 421nm (corresponding to a detuning Aq = 
— 27r X 1209GiJz from the &^P^/2 state) is optimized to 
obtain the maximal ratio between energy splittin gs o f two 
spin states and the spontaneous scattering rate [23| . 

Because of the inhomogeneity of the focused laser in- 
tensity /(r), |A£'(r)| reaches a maximum at the tar- 
get atom A and decreases dramatically at neighboring 



sites. Therefore the degeneracy of hyperfine splittings 
between different atoms is lifted. By adjusting the fo- 
cused laser intensity, the differences of the energy shifts 
5 = \^E{Q)\ - |A£'(A/2)| between target atom A and 
non-target neighboring atom B can be varied and calcu- 
lated through Eq. ([5]). Here we choose 5 — lAEu be- 
cause it balances the speed and fidelity of single spin ma- 
nipulation, leading to less spontaneously scattered pho- 
tons from the target atoms in the spin-dependent focused 
lasers. To avoid excitations of atoms to higher bands of 
the optical lattice, the rise speed of the focused laser in- 
tensity should satisfy the adiabatic condition. We use 
the adiabatic condition to estimate the ramp up time of 
the focused laser to be 35^s to give a 10"'' probability 
for excitation to higher bands. 

We then change the measurement basis by applying 
a microwave 7r/2 pulse that drives a suitable rotation to 
target atoms A (Fig. 2a). The microwave is resonant with 
the hyperfine splitting between two spin states of the tar- 
get atoms A, but has a detuning larger than 5 for non- 
target atoms B. Consider a pulse with Rabi frequency 
il{t) = exp (— Wq^^) {—tf < t < tf) and parame- 
ters ujo = U.SEr/H, Qo = ISAEr/H and tf = 5/wo. 
The pulse transfers the measurement basis of the tar- 
get atoms A in 16.9/is, while the change in the quantum 
state of non-target atoms is found to be below 3 x 10"'' by 
numerically integrating the Rabi equation that describes 
the coupling between two spin states by the microwave 
pulse. The focused lasers are adiabatically turned off af- 
ter the microwave pulse. During the step (II), the proba- 
bility for spontaneous scattering of one photon from tar- 
get atoms inside the focused laser can be roughly esti- 
mated as P « J^^ {t)dt ^ 2x 10~^, where n and tj 
represent the times when the focused laser is turned on 
and off, r is the decay rate of 6P state, d is the detuning, 
and V (t) is the potential depth of the focused laser. 

The distance between any two target atoms can be as 
short as a lattice spacing. This is because the basis trans- 
fer processes for different target atoms are preformed se- 
quentially in time, i.e., the process for an atom at site 
j2 starts after the process for an atom at site ji is ac- 
complished. In the special case that sites ji and j2 are 
spatially well separated and the final basis {l^j^) , 
at two sites are the same, the transfer process can be done 
simultaneously for two sites with one microwave pulse. 

Step (III): In this step, we transfer all atoms to the 
|P = 1) hyperfine level (Fig. 2b) to avoid stray signal in 
the detection step (IV). We apply two tt microwave pulses 
to transfer all atoms at ||) first to |P = 2, nip = 0) then 
to 1 2) = |P = l,mF = 1). Another tt microwave pulse is 
then applied to transfer all atoms at ||) to ||). The tt 
microwave pulse can be implemented within 12.5/is for a 
microwave Rabi frequency SI = 27r x 40kHz. 

Step (IV): We transfer target atoms A at jk from 
II) back to II) with the assistance of the focused lasers 
and microwave pulses (Fig. 2c), using the same atom 
manipulation procedure as that described in step (II). 
We then apply a detection laser resonant with ||) 



5 



|3) = |5^P3/2 : F = 3, m = —3) to detect the probabil- 
ity of finding all target atoms at |t) (corresponding to 
the basis state | (pj-^ ) because we transferred atoms to the 
measurement basis in step (II)) (Fig. 2d). The beam 
size of the detection laser should be large enough so that 
target atoms at different sites {jk} experience the same 
laser intensity, and scatter the same number of photons 
if they are in the |t) state. The fluorescence signal (the 
number of scattered photons) is from one of the n + 1 
quantized levels, where the l-th. level {I = 0, ...,n) corre- 
sponds to states with I sites of target atoms on state |t) 
(|(/)j^)). By repeating the whole process many times, we 
obtain the probability distribution Pi, and thus the spin 
correlation function ^ {aj^ , A; = 1, ?i} via Eq. fS]). 

The scattering photons of the fluorescence signals come 
mostly from the target atoms A at state |t) (Fig. 3). Sig- 
nal from atoms at any \F = 1) state is suppressed be- 
cause of the large hyperfine splitting {v ^ 2tt x 6.8GHz) 
between \F = 1) and \F — 2) states. Atoms B do not 
contribute to the fluorescence signal because they are al- 
ready transferred to the |F = 1) state in step (III). The 
dynamics of photon scattering is described by the optical 
Bloch equation, [23| 



dt 
d 
di 
d 
dt 



P33 



Pl3 



P31 



T 

r 



{Pl3 
+ ii^o^ Pl3 

- Wo I P31 



P3l), (10) 

^r!(2p33-l), 
^r!(2p33-l), 




10 20 30 40 50 

rt 



12 3 

n/r 



Figure 3: (a) Time evolution of the probability for the target 
atoms A to be in the excited state |3). (6) The number of 
scattering photons rip versus time for atoms A at state ||). (c) 
The number of scattering photons rip versus time for atoms 
{A or B) at states [j) and |2). (d) The number of scattering 
photons Up versus the Rabi frequency of the resonant laser for 
atoms A at state If), (e) The number of scattering photons 
rip versus the Rabi frequency of the resonant laser for atoms 
(A or B) at states |J,) and |2). The time period for (d) and 
(e) is 50/r. 



where ^33 (t) = |c3 {t)\ is the probability for the atom 
to be in the excited state |3), |1) represents state \1) for 
target atom A at state |t) and \F = 1) for other cases, 
Pi3 {t) = ci {t) cl (t) e"'"', (t) = Pis {t). The detuning 
of the laser vq is zero for target atoms A at state |t) and 
^ 1/ for all non-target atoms B and part of the target 
atoms A at hyperfine state |2). F = 27r x 6.07MHz is 
the decay rate of the excited state |3), fl is the Rabi 
frequency of the resonant laser that is related to the on- 
resonance saturation parameter by Sq — 2 |r2|^ /F^. 

We numerically integrate the optical Bloch equation 
and calculate the number of scattering photons 

np{t) = r [ P33it')dt' (11) 
Jo 

for both target and non-target atoms. In Fig. [3^, we 
plot the probability pf^ for target atoms A at state |3) 
with respect to time. We see pf^ increases initially and 
reaches the saturation value sq/2 {1 + sq). The num- 
ber of scattering photons reaches 20 in a short period, 
1.3/xs for atoms A at |t) (Fig. but it is only 10~^ 
for non-target atoms B and target atoms A at state |2) 
(Fig. El^c)). Therefore the impact of the resonant laser 
on the non-target atoms B can be neglected. In Fig. ^d) 
and (e), we see that for a wide range of Rabi frequencies 
(laser intensities), the scattering photon number for the 
non-target atoms B is suppressed to undetectable levels, 
below lO"'^. 

Unlike the noise correlation method, the accuracy of 
our detection scheme does not scale with the number of 
total atoms N, but is determined only by manipulation 
errors in the above steps. We roughly estimate that n- 
spin correlations can be probed with a cumulative error 
< n X 10^^, which is sufficient to measure both long 
and short range spin correlation functions. Our estimate 
takes into account possible experimental errors in all four 
steps, and the fact that the same experimental procedure 
is repeated many times to determine the probability Pi. 
The total failure probability of the four step detection 
process is p (n) ^ n x 10^^ for n target atoms to give 
an incorrect measurement of the target atom quantum 
state for determining the quantity Pi. Assuming that 
the experimental procedure is repeated v times, the to- 
tal probability for obtaining one incorrect measurement 
result is about vp (n) . Since one incorrect measurement 
result leads to an error ~ l/v in Pi, the expectation for 
the error is vp (n) x 1/ v — p (n) , which should be cho- 
sen to be 1/i^ to minimize the total error. In addition, 
the uncertainty in measurements of probability Pi itself 
in repeated experiments is also about l/i>. Taking into 
account of all these errors, we find that n-spin correla- 
tions can be probed with an error that scales as Cp (n) , 
where C ~ 3 in our rough estimate. As a conservative 
estimate, we take C = 10 to give an overall error in mea- 
suring n-spin correlation function ^ {aj^ , A; = 1, 71} to 
be less than lOp (n) ^ n x 10^^. 

We note that the scheme requires repeated production 
of nearly the same condensate and repeated measure- 



merits, two standard techniques which have been reaHzed 
in many experiments. We have proposed a powerful tech- 
nique for investigating strongly-correlated spin models in 
optical lattices and now consider one of several possible 
applications. 

IV. SPIN WAVE DYNAMICS: 

Our technique can be used to investigate time- 
dependence of correlation functions. In the following, we 
show how our scheme can be used to engineer and probe 
spin wave dynamics in a straightforward example, the 
Heisenberg XX model realized in optical lattices with 
a slightly different implementation scheme than the one 
discussed in the previous section. Consider a Mott in- 
sulator state with one boson per lattice site prepared in 
the state |0) = |F = 1, mp — —1) in a three dimensional 
optical lattice. By varying the trap parameters or with a 
Feshbach resonance, the interaction between atoms can 
be tuned to the hard-core limit. With large optical lattice 
depths in the y and z directions, the system becomes a 
series of one dimensional tubes with dynamics described 
by the Bose-Hubbard Hamiltonian: 

j 

This Bose-Hubbard model can be solved exactly. It offers 
a testbcd for spin wave dynamics by mapping the Hamil- 
tonian (fT2|l onto the XX spin model, ffy 2y: A = 
0) , with the Holstcin-Primakoff transformation [2a| • 

We now study the time dependent behavior of the XX 
model using our proposed scheme. In the Heisenberg pic- 
ture, the time evolution of the annihilation operator can 
be written as: aj (t) = '^j, aji (0) i-' ~^Jji-j (a), where 
Jj'-j (a) is the Bessel function of the interaction parame- 
ter a{t) = 2 Jq X {t') dt' . To observe spin wave dynamics, 
we first flip the spin at one site from | to |, which, in the 
bosonic degrees of freedom, corresponds to removing an 
atom at that site. Because of the spin-spin interactions, 
initial ferromagnetic order gives way to a re-orientation 
of spins at neighboring sites which propagates along the 
spin chain in the form of spin waves. This corresponds 
to a time dependent oscillation of atom number at each 
site. Therefore, spin wave dynamics can be studied in one 
and two point spin correlation functions by detecting the 
oscillation of the occupation probability at certain sites 
and the density-density correlator between different sites, 
respectively. 

Single atom removal at specific sites can be accom- 
plished with the assistance of focused lasers. With a com- 
bination of microwave radiation and two focused lasers 
(propagating along y and z directions respectively), we 
can selectively transfer an atom at a certain site from 
the state |0) to the state |4) = |F = 2,mF = -2). 
A laser resonant with the transition |4) |3) = 

|5^P3/2 '■ F = 3,m = —3) is then applied to remove an 




Figure 4: Site occupation probability (a) and density-density 
correlation (b) with respect to scaled spin interaction param- 
eter a oc time for A'' = 30 and: j = ((a)-solid line); j = 1 
((a)-dashed line); j = 0, j' — 3 ((&)-dashed line); and j = 2, 
j' — 3 ((fe)-solid line). The site index j ranges from — A'^ to 
iV. 

atom at that site. Following an analysis similar to the 
one above, we see that the impact on other atoms can 
be neglected. To observe fast dynamics of spin wave 
propagation, we may adiabatically ramp down the op- 
tical lattice depth (and therefore increase x) from the 
initial depth Vq = 5QEr to a final depth ISEji, with a 
hold time, ihoid, to let the spin wave propagate. Finally, 
the lattice depth is adiabatically ramped back up to Vq 
for measurement. The time dependence of the lattice 
depth in the ramping down process is chosen to be 

V (t) = Va/ (l + A^2P,^,Vo/ERLUrt^ , (13) 

where Poxc is the probability of making an excitation to 
higher bands and uju = En/h. For Poxc = 4 x lO""*, we 
find the interaction parameter to be a(ihoid) = 0.0146 -|- 
0.0228ajflihoid, with the tunneling parameter: 

X (t) = (4/V^) <V3/4 (t) exp (^-2^Vit)/En) . 

. 

Two physical quantities that can be measured in ex- 
periments are the single atom occupation probability 

D, (a) = {ip\ a+a, W) = J^, (a) (15) 

at the site j, and the density-density correlator 

Gn' (a) = ((ysl a^^ajo^^aj, |(^) (16) 
= ^ J,%.(a)J^%,(a)-^(J,%.(a)J,%.,(a) 

+Ji-j {a) J^-j (a) J^^y (a) Ji-y {a)) 

between sites j and j', where ip is the initial wavefunction 
with one removed atom at site k. The former is related 
to the local transverse magnetization through 

{^\ s] (a) 1^) = D, (a) - 1/2, (17) 

and the latter is related to the spin-spin correlator via 

G,y(a) = (^|4(a)4(a)|^) (18) 
+ {D,{a) + Dy (a))/2+l/4. 
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In Fig. 4, we plot Dj (a) and Gjj' (a) with respect to 
the interaction parameters a (which scales linearly with 
holding time). Here the initial empty site at j = is 
located at the center of the spin lattice. We see different 
oscillation behavior at different sites. Initially, all sites 
are occupied except j = 0, i.e., Dj^o (0) = 1, Dj=Q (0) = 
(Fig. 4(a)). The initial spin-spin correlation Gjj' (0) 
between different sites is zero if either j or j' is zero, 
and one if both of them are nonzero (Fig. 4(b)). As a 
increases, atoms start to tunnel between neighboring sites 
and the spin wave propagates along the one dimensional 
optical lattice, which is clearly indicated by the increase 
(decrease) of the site occupation at j = {j ^ 0). In 
Fig. 4, the oscillation of (a) and G^y (a) in a long 
time period and the decay of the oscillation amplitudes 
originate from the finite size of the spin lattice, which 
yields the reflection of the spin waves at the boundaries. 

To probe the single site occupation probability _Dj (a) , 
we use two focused lasers and one microwave pulse to 
transfer the atom at site j to |4). A laser resonant with 
the transition |4) |3) is again apphed to detect the 
probability to have an atom at |4), which is exactly the 
occupation probability Dj {a). To detect Gj^i (a), we 
transfer atoms at both sites j and j' to the state |4) and 
use the same resonant laser to detect the joint probability 
for atoms at |4). The fluorescence signal has three lev- 



els, which correspond to both atoms Gjji (a), one atom 
Dj {t) + Dji (t), and no atoms at state |4). A combi- 
nation of these measurement results gives the spin-spin 
correlator {<^\ Sj (a) s|, (a) \(p). 



V. CONCLUSION 

We find that a relation between general spin correla- 
tion functions and observable state occupation probabil- 
ities in optical lattices allows for quantitative measure- 
ments of a variety of spin correlators with the help of 
local probes, specifically focused lasers and microwave 
pulsing. Our proposal includes a realistic and practical 
quantitative analysis suggesting that n-spin correlations 
can be probed with an error < nx 10^^, which is sufficient 
to measure both long and short range spin correlation 
functions. Our work establishes a practical and workable 
method for detecting n-spin correlations for cold atoms 
in one, two or three dimensional optical lattices. Applica- 
tions to a broad class of spin physics including topological 
phases of matter [l^ realized in spin-optical lattices 
are also possible with our proposed technique. 

This work is supported by ARO-DTO, ARO-LPS, and 
LPS-NSA. 
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